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Abstract. We present optical ii-band light curves of the gravitationally lensed quasar SBS1520+530 derived from 
data obtained at the Nordic Optical Telescope. A time delay of 130 ± 3 days (la) is determined from the light 
curves. In addition, spectra of SBS1520+530 obtained at the Keck Observatory are spatially deconvolved in order 
to extract the spectrum of the faint lensing galaxy, free of any contamination by the light from the bright quasar 
images. This spectrum indicates a lens redshift z = 0.717, in agreement with one of the absorption systems found 
in the quasar spectra. The best mass model of the system includes a second nearby galaxy and a cluster of galaxies 
in addition to the main lensing galaxy. Adopting this model and an Q = 0.3, A = 0.7 cosmology, our time-delay 
measurement yields a Hubble constant of Ho = 51 ± 9km s _1 Mpc -1 (la error). 
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1. Introduction 

Prompted by the successful optical measurements of time 
delays in PG 1115+080 (Schechter et al. |1997| ), a pho- 
tometric monitoring campaign of gravitationally lensed 
quasars has been carried out at the Nordic Optical 
Telescope (NOT) with the aim of measuring time delays 
for many gravitationally lensed quasars. This observable 
is of crucial importance to study the mass distribution in 
lensing galaxies and to determine Hq using the method 
proposed by Refsdal (1964). 

Two time delays have been measured so far at the 
NOT: one for the doubly imaged quasar B 1600+434, 
lensed by an edge-on spiral (Burud et al. 2000) and an- 
other for the two (summed) components RX J0911+0550 
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Space Administration. 



(Hjorth et al. 2002). In this paper we will present our 
third time delay measurement: that of SBS1520+530, a 
doubly imaged BAL quasar at z = 1.86. SBS 1520+530 



was discovered by Chavushyan et al. ( |1997| ) as a dou- 
ble quasar with angular separation of 1.56". The lensing 
galaxy was detected by Crampton et al.(1998) on ii-band 
images obtained by using the adaptive optics system of 
the Canada-France-Hawaii Telescope. The lensing galaxy 
is also well resolved on optical images obtained at the NOT 
and deconvolved using the MCS deconvolution algorithm 
(Magain, Courbin & Sohy 1998), and on public Hubble 
Space Telescope near-IR data (Faure et al. 2002). The 
redshift of the galaxy however, remains unknown. With 
the aim of measuring this redshift, we have obtained a 
spectrum of SBS 1520+530 with ESI at the Keck obser- 
vatory. 

The observations and data reduction of the images are 
presented in Sect. ^| and |3|. The time delay measurement 
is described in Sect. |] while the analysis of the spectro- 
scopic data is explained in Sect. [g. This section also in- 
cludes a discussion on the spectral differences between the 
two quasar images. Mass models and the estimate of the 
Hubble constant are discussed in Sect. |. Finally, Sect. @ 
summarises the main results. 
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Fig. 1. Deconvolved i?-band image (FWHM=0'/28) of 
12"x 12"around SBS1520+530 obtained from the simulta- 
neous deconvolution of 59 frames. We see the two stars just 
to the North- West (S2) and South-East (SI), the quasar 
components A and B, the lensing galaxy L close to B and 
an additional galaxy M, NE of the lens. North is up and 
East is to the left. 

2. Observations and Data Reduction 

Weekly observations of SBS 1520+530 were carried out 
at the NOT from February 1999 to May 2001. The de- 
tector ALFOSC (AndaLucia Faint Object Spectrograph 
and Camera) with a pixel size of 0'.'188 was used for 95% 
of the frames obtained. The remaining 5% were obtained 
with HiRAC (High Resolution Adaptive Camera) or with 
the stand-by camera StanCam which have pixel sizes of 
0'.'107 and 0'.'176 respectively. 

The target is relatively bright (I~18) and an exposure 
time of 300 sec in the i?-band was sufficient to obtain a 
S/N^100 for each quasar image. The seeing varied from 
0'.'5 to 2'.' 2 during the two years of observations, l'.'O being 
the most frequent value. 

Our automated pipeline employing routines in the 
IRAF/NOAO/CCDRED package was used in order to 
pre-process the CCD frames in an efficient and homoge- 
neous way. 



3. Photometry 

The photometry of the blended quasar images was per- 
formed by applying the MCS deconvolution algorithm 
(Magain, Courbin & Sohy 1998 ). This algorithm has al- 



Fig. 2. The field (~5'x5') around SBS1520+530 observed 
with ALFOSC at the NOT . North is up and East to the 
left. See Faure et al. 2002 for a detailed investigation of 
the spatial distribution of secondary lenses in this field 



all the frames from different epochs simultaneously, hence 
constraining the astrometry of the two quasar images and 
the shape of the lensing galaxy using the total S /N of the 
whole data-set. The intensity of the point sources are al- 
lowed to vary from image to image, which produces the 
light curves. 

The two quasar components are clearly distinguished 
on our deconvolved image of SBS1520+530 (see Fig. 0), 
and the lensing galaxy is resolved. The lensing galaxy can 
easily be taken into account in order to avoid contamina- 
tion of the photometry of the quasar images. 

The light curves of SBS1520+530 consist of 57 data 
points in the i?-band as presented in Fig. || and Table |[ 
In addition to the two stars just to the South-East (SI) 
and North- West (S2) of the quasar (Fig. |l]), three other 
reference stars in the field were deconvolved, in order to 
check the relative photometry and to check for system- 
atic errors (stars 3 to 5 in Fig. ||). The magnitudes of the 
quasar are relative to the magnitudes of Star 3 on Fig. |[ 
The error-bars include both photon noise and additional 
systematic errors, e.g., PSF errors. The latter is estimated 
by using the reference stars, as explained in Burud et al. 

( poool ). 



ready been successfully applied to several moni toring data 
of lense d quasars (e.g., Burud et al. 2000 , Burud et 
al. |2002| ). Its main advantage is its ability to deconvolve 



4. Time delay measurement 

By sliding the light curves across one another one can 
derive a rough time delay estimate of ~ 125 days. As pre- 
dicted by gravitational lens theory, the brightest quasar 
image A is the leading component. Using the \ 2 mini- 
mization method described in Burud et al. a time 
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Fig. 3. i?-band lightcurves of SBS1520+530 and of the 
reference star S4 (see Fig. ||), calculated relative to the 
reference star S3 (Fig. §. The error bars include photon 
noise and PSF errors estimated from the deconvolution of 
a reference star. The B component is shifted vertically by 
—0.4 mag and star S4 is shifted by +0.5 mag on the plot. 



delay value of 128 ± 3 (la) days is found from the i?-band 
light curves. The magnitude difference between the A and 
the B component is 0.692 ± 0.003 mag, corresponding to 
an A/B flux ratio of 1.892 ± 0.002. The error bar on the 
time delay is obtained from Monte Carlo simulations of 
1000 sets of light curves, assuming that the photomet- 
ric errors are uncorrelated and follow a Gaussian distri- 
bution. The best x 2 -fit is obtained when a linear term, 
modelling slow microlcnsing effects, is added to one of 
the components. The slope of this term is determined by 
the algorithm. With this external variation, the shifted 
light curve of the B component plotted in Fig. || becomes 
B' = B - 0.692 - 0.000067 x (JD - 2451226). The lin- 
ear correction does not remove all the external variations. 
Faster variations on time scales of 50 — 100 days are still 
present. The iterative version of the algorithm (Burud et 
al. " 



2001) was also applied in an attempt to correct for 



these fast variations. This method yields At = 130 ± 3 
days, slightly higher than the value found with the direct 
method but still in agreement within the error bars. Even 
with the iterative method some of the fastest variations 
are not corrected for (Fig. notably the small peak in 
the A component at — JD 2451750 (abscissa 850 in Fig. 0). 



5. Spectroscopy of SBS1520+530 



Fig. 5. Spectrum of the lensing galaxy in SBS1520+530. 
The Call (A 3933, 3968 A) is identified. Also indicated are 
a night sky emission line and an atmospheric absorption 
line. 



5.1. Keck II spectroscopy 

Converting the time delay into an estimate of the Hubble 
parameter requires good knowledge of the geometry of 
the lensing system, including the redshift of the source 
and of the lensing galaxy. For this purpose, we obtained 
a deep optical spectrum of SBS1520+530 with the Keck 
II Telescope and th e Ech ellette Spectrograph and Imager 
(ESI) (Sheinis et al. 2000 ). The observations were obtained 
on May 23, 2001 under good seeing conditions (0.6"). The 
1" slit was aligned to obtain simultaneously the spectrum 
of SBS1520+530 and of star S2 (indicated on Fig. 0). Four 
900s exposures were taken so that cosmic rays could be ef- 
ficiently removed. 



5.2. Redshift of the lensing galaxy 

Using the spatial information contained in the spectrum 
of the PSF star S2, the 2D spectrum of SBS 1520+530 
was spatially deconvolved, by applying the spectral ver- 
sion of the MCS deconvolution algorithm (Courbin et 
al., 2000b). The deconvolution process decomposes the 
data into the individual spectra of the two quasar images 
and of the faint lensing galaxy (see Figs. || and ||), as 
done for the lensed systems HE 2149—2745 (Burud et al. 
|2002D or HE 1104-1805 (Lidman et al. |2000| ). No emis- 
sion lines were detected in the spectrum of the lensing 
galaxy (Fig. ||). However, an absorption doublet is clearly 
observed at A ~ 6700 A, matching within the error-bars 
the Call H and K lines at redshift 0.71. This strongly sup- 
ports a scenario where the lensing galaxy is responsible for 
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Fig. 4. Time delay shifted light curves of SBS1520+530. Left: B is shifted by 0.69 mag, no correction for external 
variations (e.g., as can be induced by microlensing) is applied. Middle: B is shifted by 0.69 mag and a linear correction 
for the external variations is applied. Right: B is shifted vertically by 0.69 mag and corrected for external variations 
with the iterative method (cf. Sect. 0). 



the CaH absorption lines at redshift z = 0.71 ± 0.005 first 



found by Chavushyan et al. (1997) in the quasar spectra. 
Another absorption system was found by Chavushyan et 
al. (1997), at z = 0.82. This system does not seem to be 
associated with the lens. 



5.3. Spectra of the quasar components 

The quasar spectra display the usual prominent emission 
lines such as A1III/CIII] (A 1909 A) and Mgll (A 2798 A), 
at redshift z = 1.855. Both (deconvolved) spectra show 
the metal absorption lines first reported by Chavushyan 
et al. (1997). These are the Fell (A 2587, 2600 A), Mgll 
(A 2798A), Mgl (A2853 A) absorption lines at redshifts 
0.717 ± 0.001 and 0.817 ± 0.001. Call (A 3933, 3968 A) 
is seen only at z = 0.717 and is indicated on Fig. |^. It is 
better seen in the normalized spectra of Fig. ^| 

The general shape of the spectra towards the blue 
wavelengths is most probably corrupted by atmospheric 
refraction. The observed reddening of B compared to A 
at the blue wavelengths may therefore not be real. In the 
red parts of the spectra, we find a flux ratio of 1.82 (at A 
9000 A) which, compared with the ratio of 1.892 measured 
on R band light curves (A 6500 A), indicates an almost 
constant flux ratio in this wavelength range. 

The continuum normalized spectra, displayed in Fig. ^, 
show that the equivalent widths of the emission lines are 
much larger in A than in B, suggesting that B's continuum 
is amplified. Such a differential amplification can be ex- 
plained by microlensing effects. Assuming that the lensed 
quasar follows the "standard" AGN model and that we 
are observing component B through a network of caus- 
tics produced by stars in the main lensing galaxy, one can 
imagine that the inner continuum region of component B 
is enhanced by a larger amount than the much larger emis- 
sion line region. This interpretation has already been pro- 
posed to explain similar spectral differences observed in 



the d ouble H E 1104-1805 (Wisotzki et al. |1993|, Co urbin 
et al. |2000bD and HE 2149-2745 (Burud et al. |2002fi. 



6. Modelling the lens 

The time delay measured for SBS 1520+530 can be used 
to infer an estimate of Hq, via a model of the total grav- 
itational potential responsible for the lensing effect. This 
model shall take into account the main lensing galaxy and 
a small galaxy cluster located to the North- West of the 
quasar. 

Our models, already described in Faure et al. (2002), 



assume that different mass components along the line of 
sight are Pseudo-Isothermal-E lliptic al-Mass-Distri bution 
(PIEMD, K assiol a & Kovner |l993[ Kneib et al. |l993 , 
Kneib et al. 1996). In addition to the obs ervat ional con- 
straints already available in Faure et al. ( |2002 ), we now 



also try to pin down the flux ratio between the two quasar 
components, and to use it in the models. This ratio is dif- 
ficult to estimate for SBS 1520+530: it varies with wave- 
length, it has different values in the continuum and in the 
emission lines, and even varies within the emission lines. 
Microlensing shall be blamed, but one should also keep in 
mind that our spectra were not taken at parallactic angle, 
with a relatively narrow slit, and at large airmass. The 
values we derive are therefore only tentative: tub - niyl = 
1.40 mag (i.e., Fa/Fb — 3.6) in the emission lines, and 
tl\b - niyi = 0.83 mag (i.e., Fa/Fb — 2.15) in the contin- 
uum The latter value was measured from X-band images 
(Faure et al. |2002| ). 

Table [l] presents two lens models, computed for the 
two flux ratios, and matching a time-delay of 130+3 days. 
Model (L) considers only the main lensing galaxy L (see 
[l]) with a measured ellipticity e = 0.50 ± 0.2 and with 
a position angle PA = — 23deg±7deg (see Faure et al. 
2002). Model (L+M+C) includes the main lensing galaxy 
(L), galaxy M, located at 2.6" North-East from the quasar 
image A (see Fig. Q), and a galaxy cluster, centered at 
~1' to the North- West of the quasar. 
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Table 1. Lens models: the first column describes the lens parameters, where 7 and 9 1 are the intensity and orientation 
of the shear at the position of the quasar images. Column (2) and (3) gives the models results for the two possible 
flux ratio measured in SBS f 520+530. The letter a refers to lens velocity dispersions and r CMt is the cut radius of the 
different mass distributions. 



Mb-Ma 


0.83 mag 


1.4 mag 


L 


no good fit 


az,=228 km s _1 r cll t(L)=13 kpc 


7 




0.39 
+75.3 deg 


Ho (la) 




63 ± 9 km s^Mpc" 1 


L+M+C 


ai, = 189 km s" 1 r cui (L)=51 kpc 
a M = 135 km s _1 r cut (M)=20 kpc 
ac=718 km s _1 r cu t(C)=600 kpc 


a L =200 km s" 1 r cut (L) = 100 kpc 
a A /=118 km s _1 r cut (M)=40 kpc 
a c =400 km s _1 r cut (C)=600 kpc 


7 


0.32 
+74.5 deg 


0.34 
+74.3 deg 


Ho (la) 


52 + 8 km s-'Mpc -1 


50 + 8 km s-'Mpc -1 



Modeling SBS 1520+530 is possible with a single lens- 
ing galaxy (i.e., galaxy L), only if - m^ > 1.15. Given 
the errors on the image flux ratio, such a model is plausi- 
ble. Taking the emission line flux ratio (column 3 in Table 
|l|) as the true one, because it is less affected by microlens- 
ing, we find Ho = 63±9 km s" 1 Mpc^ 1 . The error, once the 
image flux ratio is fixed, comes mainly from the position 
angle of the lens and its ellipticity. 

If the magnitude difference between the images is less 
than 1.15, SBS 1520+530 can not be modeled by usin 
one single galaxy. The model in column (2) of Table [l| 
with galaxy M and an intervening galaxy cluster gives 52 
± 8 km s _1 Mpc" 1 . 

Finally, changing the flux ratio does not change much 
the value of Ho, for a given lens model. The flux ratio is 
important in the choice of model (i.e., multiple or single 
lens) but does not alter much the fit once the model is 
given. 



7. Discussion 

The time delay At — 130±3 days (1 a) has been measured 
for the first time, in the lensed quasar SBS 1520+530 on 
the basis of i?-band images obtained with the NOT. Keck 
spectroscopy of the lensing galaxy strongly suggests that 
the absorption system at z — 0.717 is associated with the 
lensing galaxy. 

Applying t he d etailed mass model presented by 
Faure et al. (2002) we derive a mean value Ho = 
51 ± 9 kms _1 Mpc _1 . When only the main lensing 
galaxy is used in the mass model Ho increases to 63 
±9 kms _1 Mpc— 1. The fit of the wo different models de- 



pends on the flux ratio between the two images, the model 
with only one galaxy only gives a good fit when the high 
flux ratio measured from the emission lines is applied. 
Determining precisely the flux ratio is therefore crucial to 
further improve the modeling of the system and hence the 
precision on Ho. Finally, a possible systematic error adds 



to the uncertainty in Ho if the central mass concentration 
is not isothermal. 

External variations, probably due to microlensing ef- 
fects, are observed on the time delay shifted light curves 
(Fig. ||). Part of this effect can be modeled as a linear term, 
or corrected for with the iterative algorithm for measur- 
ing time delays. However, significant external variations of 
time scales of ^50 days do remain once these corrections 
are made. 

Microlensing is important in SBS 1520+530 as sug- 
gested both by the high frequency variations in the light 
curves and by the spectra of the quasar images. With light 
curves in only one band we can not efficiently disentan- 
gle between microlensing and intrinsic variations of the 
quasar. Such a work shall be possible using the colour in- 
formation provide by a monitoring of the object in spec- 
troscopy. The issue in conducting such a spectrophotomet- 
ric monitoring for SBS 1520+530 is double: (1) to mea- 
sure the true flux ratio between the quasar images and dis- 
criminate between lens models and (2) to use microlensing 
to infer constraints on the distribution of micro-lenses in 
galaxy (L) and/or to reconstruct the energy profile of the 
central AGN in the source. 
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Fig. 7. The A (top) and B (bottom, vertically shifted 
by -0.2 for clarity) spectra with normalized continuum. 
The equivalent widths of the emission lines including 
A1III/CIII] and Mgll are significantly smaller in B than in 
A. The CaH absorption lines at redshift 0.7f are indicated 
as well as two telluric bands. 
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Table 2. R-band photometry of SBS 1520+530 



A 






B 




JD 


1 


.2548 


.0143 


1. 


.868 


.0242 


2451226.75 


1 


.1809 


.0091 


1. 


.865 


.0115 


2451235.75 


1 


.2221 


.0087 


1. 


.848 


.0096 


2451243.75 


1 


.2245 


.0087 


1. 


.893 


.0097 


2451255.75 


1 


.1905 


.0086 


1. 


.893 


.0096 


2451264.75 


1 


.2312 


.0086 


1. 


.935 


.0095 


2451284.50 


1 


.2257 


.0087 


1 


.943 


.0098 


2451306.50 


1 


.2154 


.0087 


1 


.937 


.0101 


2451311.50 


1 


.1879 


.0087 


1 


.911 


.0097 


2451319.50 


1 


.1747 


.0086 


1 


.931 


.0096 


2451336.50 


1 


.1725 


.0087 


1. 


.922 


.0099 


2451342.50 


1 


.1629 


.0087 


1 


.915 


.0097 


2451354.50 


1 


.1426 


.0085 


1. 


.905 


.0091 


2451362.50 


1 


.1448 


.0089 


1. 


.915 


.0109 


2451371.50 


1 


.1341 


.0086 


1. 


.916 


.0092 


2451377.50 


1 


.1426 


.0086 


1. 


.921 


.0096 


2451386.50 


1 


.1525 


.0087 


1 


.931 


.0102 


2451403.50 


1 


.1515 


.0086 


1. 


.963 


.0095 


2451420.25 


1 


.1293 


.0085 


1. 


.922 


.0090 


2451433.25 


1 


.1310 


.0085 





.920 


.0089 


2451434.25 


1 


.1341 


.0089 


1. 


.880 


.0096 


2451443.25 


1 


.1406 


.0086 


1. 


.899 


.0091 


2451450.25 


1 


.1449 


.0087 


1 


.865 


.0099 


2451457.25 


1 


.1433 


.0089 


1 


.850 


.0107 


2451465.25 


1 


.0952 


.0087 


1 


.833 


.0096 


2451530.75 


1 


.1086 


.0086 


1. 


.822 


.0094 


2451566.75 


1 


.1139 


.0087 
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